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Abstract 

We argue that double logarithmic corrections a s In 2 x need to be resumed in perturbative QCD factorization theorem 
for exclusive B meson decays, when the end-point region with a momentum fraction x — > is important. These double 
£SJ 1 logarithms, being of the collinear origin, are absorbed into a quark jet function, which is defined by a matrix element 
of a quark field attached by a Wilson line. The factorization of the jet function from the decay B — > jW is proved to 
all orders. Threshold resummation for the jet function leads to a universal, i.e., process- independent, Sudakov factor, 
whose qualitative behavior is analyzed and found to smear the end-point singularities in heavy-to-light transition form 
factors. 



o 



I. INTRODUCTION 



Perturbative QCD (PQCD) factorization theorem for exclusive processes |T) states that a hadronic form factor 
at large momentum transfer is expressed as convolution of a hard amplitude with hadron distribution amplitudes. 
However, the PQCD evaluation of the pion form factor suffers the soft enhancement from the end point of a momentum 
fraction x — > |2| . In this region the hard amplitude is characterized by a low scale, such that perturbative expansion 
in terms of a large coupling constant a s is not self-consistent. More serious end-point singularities, logarithmic and 
linear, have been observed in the leading-twist (twist-2) and next-to-leading-twist (twist-3) contributions to the B — ► 7r 
t— I i transition form factors ||-||], respectively. 

We argue that when the end-point region is important, the double logarithms a s In 2 x from radiative corrections 
should be organized to all orders in order to improve perturbative expansion. These double logarithms have been 
O |' found in the semileptonic decay B — > irlD Q| and in the radiative decay B — > jlD jfj), and resummed. Here we shall give 
a systematic treatment of this type of double logarithms, which is of the collinear origin, by introducing a quark jet 
function into PQCD factorization theorem. The all-order factorization of the jet function from the decay B — ► jW is 
proved following the procedures proposed in [[7||| , which provides a solid theoretical ground for the modified formalism 
appropriate for the end-point region. It will be shown that the jet function is defined as a matrix element of a quark 
. ' field attached by a Wilson line. 

The double logarithms in the et function are resummed into a Sudakov form factor in the Mellin space using 
' the technique developed in [fj|-[l2|. In principle, with the definition of the jet function constructed in this work, 
its threshold resummation can be performed up to the next-to-leading-logarithm accuracy. Since our purpose is to 
understand the qualitative behavior of the jet function at small x, it suffices to derive only the leading-logarithm result. 
Extending the above formalism to the semileptonic decay B — > irW, we obtain the identical Sudakov factor, implying 
its universality. By means of the inverse Mellin transformation, the Sudakov factor is found to vanish quickly as x — > 0. 
The suppression is so strong that the end-point singularities in the B — > ir form factors are smeared. We conclude 
that in a self-consistent perturbative evaluation of heavy-to-light transition form factors, the end-point singularities 
do not exist. 

It has been argued that the dependence on the transverse momentum kx carried by the light spectator in a B meson 
can not be dropped from hard amplitudes for exclusive B meson decays ||. As considering higher-order corrections 
to a hard amplitude, kx appears in the ratio kx/k + , k + being the longitudinal component of the light spectator 
momentum. Obviously, this ratio is not power-like, because kx and k + are of the same order of magnitude. Hence, 
a reliable formalism for exclusive B meson decays must involve the parton transverse degrees of freedom. The effect 
due to the inclusion of kx has been explored thoroughly in fil3-[lB], which also smears the end-point singularities in 
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the heavy-to-light form factors [[L6|. In this paper we shall concentrate only on the Sudakov effect from threshold 
resummation, and neglect the kx dependence. 

In Sec. II we work out the 0(a s ) factorization formula for the B — > decay amplitude at the end point. The 
all-order factorization is proved in Sec. Ill, leading to the convolution of the hard amplitude with the B meson 
distribution amplitudes and with the jet function. Threshold resummation is done in Sec. IV by solving an evolution 
equation for the jet function. In Sec. V we extract the jet function from the semileptonic decay B — > irW. The 
smearing effect from the Sudakov factor on the end-point singularity is demonstrated numerically by evaluating the 
leading-twist contribution to the B — > ir form factor. Section VI is the conclusion. In Appendix A we comment on 
the Sudakov resummation performed in H . 



II. O(as) FACTORIZATION 



We study the radiative decay B — > 7Z1/. The B meson carries the momentum P\ = (Afs/v / 2)(l, I, Or) and the 
outgoing photon carries the momentum P2 = (Mb/ V%)(0, f], Or) with the energy fraction rj, where the light-cone 
variables have been adopted. Consider the kinematic region with small q 2 = (1 — rf)M B , q = P\ — P2 being the lepton 
pair momentum, i.e., with large 77, where PQCD is applicable. The mass difference between the B meson and the 
b quark, A = Mg — nib, is treated as a small scale in the heavy quark limit. As stated in [Q, there are two types 
of infrared divergences in radiative corrections, soft and collinear. Soft divergences come from the region of a loop 
momentum I, where all its components vanish. Collinear divergences are associated with a massless parton carrying 
a momentum of order Mb- The 0(a s ) factorization of the soft divergences into the B meson distribution amplitudes 
has been presented in 0. Here we derive only the 0(a s ) factorization of the collinear divergences at the end point 
into the jet function. 

According to leading-twist PQCD factorization theorem, the form factor relevent to the decay B — ► is written 

as 



(1) 



where the symble ® denotes the convolution over the spectator momentum fraction x = k + /P^, k being the spectator 
quark momentum. The light-cone B meson distribution amplitudes <f) m are defined by 



\v(y )75(/' + ^)7 ± exp 



dzri- ■ A(znJ) 



b v {Q)\B{Pi)) , 



(2) 



with the dimensionless vectors v — P\/Mb and n_ = (0, 1, Ot), and the rescaled b quark field b v . The hard amplitudes 
H m are obtained by contracting the quark- level diagrams with the spin structures corresponding to 4> m [ pT| . Figure 
1 is the lowest-order example, where the upper line represents a b quark and x represents a weak decay vertex. The 
contributions to the hard amplitudes from Figs. 1(a) and 1(b) scale like 1/(AMb) and 1/M B , respectively. Below we 
shall concentrate on Fig. 1(a). 

The factorization formula in Eq. ([l]) is appropriate for the region with k + ~ 0(A), in which the only infrared 
divergences are the soft ones absorbed into <f>± in Eq. (g). Near the end point k + ~ 0(A 2 /Mb), the internal quark in 
Fig. 1(a) carries a large momentum P^ — k with its invariant mass vanishing like (P2 — k) 2 — —2xP\ ■ P2 ~ (9(A 2 ). This 
kinematics is similar to the threshold region of deeply inelastic scattering (DIS) with the Bjorken variable Xb — * 1, 
where the scattered quark also carries a momentum with a large component and possesses a small invariant mass 
(1 — xb)s, s being the center-of-mass energy. In this region the scattered quark produces a jet of particles, to which 
the radiative corrections contain additional collinear divergences. Hence, a jet function needs to be introduced into 
the standard factorization theorem for DIS |18[ | . Similarly, a jet function has been incorporated into the factorization 
of direct photon production at a large photon transverse momentum (threshold) jl9| . Here we associate a jet function 
with the internal quark at the end point of the momentum fraction involved in the decay B — > jW. 

An additional collinear divergence from the loop momentum parallel to P2 appears in the higher-order correction 
to the weak decay vertex shown in Fig. 1(c). This divergence can be extracted by replacing the b quark line by an 
eikonal line in the direction of n + : 



/Pi- jk+ / + m b 
(Pi-k + l) 2 - m\ 



7 /3 &(Pi - k) 



»+ 



b(Pi - k) 



(3) 
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where b(P\ — k) represents the b quark spinor with the momentum Pi — k and n+ = (1,0, Ot) is a dimensionless 
vector. The above replacement makes sense, since 7^ is dominated by its plus component 7 + , l~ is much larger than 
l + , and py" diminishes. The factorization of the fermion flow is achieved by inserting the Fierz identity, 



I in I] 



ij Ilk 



\hkhj + ^(75)^(75)0- + ^(7a) 4 fc(7 Q )ij + ^(757ak-(7 Q 75)ij + ^{(J a i3)ik(<J al3 )ij 



in which the first and last terms contribute in the combined structure, 



Assigning the identity matrix / in Eq. 
jh+ jti-jA. then leads to the loop integral, 



lijhk — * -^iik(A+ A-)ij 



(4) 



(5) 



to the trace for the hard amplitude, we obtain Fig. 1(a). The matrix 



-ig 2 C F 



dH 1 



(2tt) 4 4 



fi+ A-7/3 



(P 2 - k + 1) 2 



■IP ' 



4tt 



Cp In x 



(6) 



where Cp = 4/3 is a color factor. To obtain the above double logarithm, we have regularized the collinear pole from 
I parallel to n + by allowing n + to possess a small amount of minus component. The correction to the photon vertex 
in Fig. 1(d) contains only the single logarithm a s lnx, since the phase space of the loop momentum is restricted to 
< l + < k + ~ 0(A 2 /Mb)- The self-energy correction in Fig. 1(e), generating also the single logarithm, is factorized 



into J'l trivially by applying the Fierz transformation in Eq. (|5|). The explicit expression of J^' is not essential here. 
The reason for labeling the 0(a s ) jet functions by || and _L will become clear in the next section. To group Figs. 1(c) 
and 1(e), we introduce a quark jet function with the zeroth-order expression and the 0(a s ) expression, 



(i) 



= tr{jh+ >»_)/4 = 1 , J (1) = jj 13 + Jj 1 



(J) 



respectively. 

It has been shown that the soft divergences in Figs. 1(c), 1(d) and 1(f)- 1(h) are absorbed into the 0(a s ) B meson 
distribution amplitudes [Q. The remaining infrared finite 0{a s ) contributions, including the single logarithms 

r(i) 



except that from Fig. 1(e), are assigned into the hard amplitudes H±' 
end-point region is written, up to 0(a s ), as 



Therefore, the modified factorization in the 



G (o) +G (D = £ [1 + ^) ]0[if (o )+jff (i) ]0[1 + J (i) ] + o(a 2 ) 



(8) 



with G^ ) and GK 1 ' displayed in Fig. 1. Note that J does not depend on the subscript m, a fact obvious from the 
derivation of Eq. (0) . 



III. ALL-ORDER FACTORIZATION 



In this section we prove the factorization theorem for the radiative decay B — > jlv at the end point to all orders, 
and construct the definition of the jet function J(x) via 



J(x)u(P 2 - k) = (u{P 2 - fc)|g(0)i A+ A- exp 



dzn+ ■ A(zn + ) 



|0> 



(9) 



The spinor u(P 2 — k) is associated with the internal quark, through which the momeutm P 2 — k flows. The idea of 
the proof is based on induction J7|,p| with the help of the Ward identity, 



l^{l,kxM,---,k p ) = 



(10) 



where G M is a physical amplitude with an external gluon carrying the momentum I and with p external quarks carrying 
the momenta hi, k 2 , • • •, k p . All these external particles are on mass shell. 
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The 0(a s ) factorization of the collinear divergences has been worked out in Sec. II. Assume that the factorization 
theorem holds up to 0(a^): 

G= <l> m ®H m ®J, (11) 

m=+,— 

with the definitions, 

N N N N 

i=0 i=0 i=0 i=0 

G^ denotes the full diagrams of 0{a\). The B meson distribution amplitudes d>m and the jet function J^> are 
defined by the 0{a\) terms in the perturbative expansions of Eq. (0) and of Eq. (H), respectively. The 0(a s ) hard 

(i) ... . ... 

amplitudes Hm do not contain the intial-state soft divergences and the end-point double logarithms, which have been 
collected into the B meson distribution amplitudes and into the jet function, respectively. We then have the relations, 

k k — i 

G ( fe)= £ £X>« ® ®J«>, k = 0,l,---N. (13) 

m= + ,- i=0 j=0 

Below we prove the collinear factorization of the 0(af +1 ) diagrams G' JV+1 ' into the convolution of the O(a^) 
diagrams G^ N ' with the 0(a s ) jet function jW. Look for the gluon in a complete set of 0(a^ +1 ) diagrams 

G (N+1) 

one of whose ends attaches the lower most vertex on the internal quark line. Such a gluon exists, because G^" 1 " 1 ) are 
the finite-order diagrams. Let a denote the lower most vertex, and (3 denote the attachments of the other end of the 
identified gluon inside the rest of the diagrams. There are two types of collinear configurations associated with this 
gluon, depending on whether the vertex [3 is located on an internal line with a momentum along P%. The fermion 
propagator adjacent to the vertex a is proportional to pi in the collinear region with the loop momentum I parallel 
to Pi. If P is not located on a collinear line along Pz, the component 7~ in 7" and the plus component of the vertex 
(3 give the leading contribution. If j3 is located on a collinear line along P2, f3 can not be plus, and both a and (3 label 
the transverse components. This configuration is the same as of the self-energy correction to an on-shell particle. 

According to the above classification, we decompose the tensor g a p appearing in the propagator of the identified 
gluon into 

n +oth x c 1 ( n+Jp \ 

g a /3 = r - O a X_(>BX + \9al3 7 + <WtV3_L • (14) 

n+ • / \ n + ■ I J 

The first term on the right-hand side of Eq. (|lj) extracts the first type of collinear divergences, since the light-like 
vector n+ Q selects the minus component of 7°, and lp in the collinear region selects the plus component of the 
vertex (3 . The second term extracts the second type of collinear divergences. The last term does not contribute a 
collinear divergence. We shall concentrate on the factorization corresponding to the first term, and the factorization 
corresponding to the second term can be achieved simply by requiring gauge invariance |s| . 

The identified collinear gluon with a = — and (3 — + does not attach the internal quark line directly, which carries 
a momentum along Pi at the end point. That is, those diagrams with Fig. 1(e) as the 0(a s ) subdiagram are excluded 
from the set of G<- N+r > as discussing the first type of collinear configurations. Applying the Fierz transformation in 
Eq. (^) to break the fermion flow, we have the physical amplitude, in which the two on-shell quarks and the on-shell 
gluon carry the momenta £P2, £,P2 — I and Z, respectively. The fraction £ reflects that the momentum flowing through 
the internal quark is almost parallel to P^. Figure 2(a), describing the Ward identity, contains a complete set of 
contractions of lp represented by arrows, since the second diagram has been added back. The cuts on the internal 
quark lines denote the insertion of the Fierz identity. 

The second diagram in Fig. 2(a) gives 

where the factor £ /Pi at the beginning of the above expression comes from the internal quark propagator adjacent to 
the photon vertex, and the term — £ /Pi at the end leads to the 0(a^) diagrams. The diagrams G^ corresponding 
to the first term in Eq. ( |T^ ) are then factorized according to Fig. 2(b). The factor n +a /n + ■ I is just the Feynman 
rule associated with the Wilson line along the vector n+ in Eq. (|^). We arrive at the convolution for the first type of 
collinear configurations, 
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G (N+i) M G[N ) 8 jji) ^ (16) 

with J|j 1 ' > given in Eq. (^|). 

The above procedures are applicable to the 0(a{ +1 ) jet function J^ +1 \ We identify the gluon in a complete set 
of 0(a J s +1 ) diagrams J^ +1 \ one of whose ends attaches the lower most vertex a on the internal quark line. The 
other end attaches the vertex (3 inside the rest of the diagrams. For the first term on the right-hand side of Eq. (fill), 
there exists a Ward identity similar to Fig. 2(a). It is then trivial to obtain the factorization of the jet function 
corresponding to the first type of collinear configurations, 

J« +1 W«®J«. (17) 

The rest part of the proof is exactly the same as of the factorization theorem for the semileptonic decay B — ► ttW in 
Sec. IV B of [0, with the pion distribution amplitudes replaced by the jet function, both of which are of the collinear 
origin. Employing Eqs. (|l3|), ( |l6| ) and (|lj]) and following the steps of factorizing the decay B — > ttIv, we derive 

N+l N+l-i 

G{N+ i) = j2 E 4® « h^-d ® jw , (is) 

m=+,- i=0 i=0 

with the 0(a^ +1 ) hard amplitude Hm +1 ^ being infrared finite. Equation ( |l8| ) indicates that all the soft and collinear 
divergences at the end point in the radiative decay B — > jW can be factorized into the B meson distribution amplitudes 
and the jet function, respectively, order by order. We complete the proof of the factorization theorem, which is 
graphically described in Fig. 2(c), and construct the jet function defined in Eq. (|]). 



IV. JET FUNCTION 

We now perform threshold resummation of the double logarithms a s In 2 x in covariant gauge d • A = 0, which 
have been collected into the jet function to all orders in the previous section. Threshold resummation for inclusive 
QCD processes has been studied intensively [p|,F0L Here we shall adopt the framework developed in fL2|[ , which has 
been shown to lead to the same results as in |J10|. First, we replace the vector n + by n, which contains a plus 
component n + > and a (small) minus component. This replacement, regularizing the collinear pole, extracts the 
double logarithm as stated in Sec. II. The jet functions in terms of the dimensionless vectors n + and n can be regarded 
as being defined in the different factorization schemes. The definition then involves three variable vectors: the Wilson 
line direction n, the large momentum P%, and the spectator momentum k. We argue that the jet function depends on 
the Lorentz invariants, n ■ k, n ■ P2 and P2 ■ k, since the other invariants P 2 2 an d ^ 2 vanish, and n 2 will be fixed below. 
The third invariant P2 ■ k is not independent, because it can be rewritten as n ■ P^n ■ k with n 2 being a constant. We 
further argue that the scale invariance in n, as indicated by the Feynman rule associated with the eikonal line along 
n, implies that the jet function must depend on k through the ratio n ■ k/n ■ P2. 

The next step is to derive the evolution of the jet function in x, i.e., in k + = xP^ by considering the derivative, 

+ dJ _ n-k dJ 

k dk^-J^k P2 d^ 7 (19) 

where we have applied the chain rule to relate the derivatives with respect to k and to n. The differentiation d/dn a 
operates on the eikonal line along n, giving 

n ■ k „„ d n., h,. 
P2 ■ k dn a n ■ I n ■ I 

with the special vertex, 

n-k P 2 -l 

n » = -p2-k— n »- (21) 

Another term in Eq. (pl|), proportional to P2, is negligible, as the special vertex attaches inside of the jet function, 
which is dominated by momenta parallel to P 2 . 

The loop momentum / flowing through the special vertex does not generate a collinear divergence due to vanishing 
of the numerator P2 ■ I in this region. It is easy to confirm that the ultraviolet region of I docs not produce In a: 
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either. Therefore, we concentrate on the factorization of the soft gluon emitted from the special vertex, which can be 
achieved by applying the eikonal approximation to internal quark propagators, leading to n_„/n_ ■ I. Following the 
reasoning in pa, the derivative of the jet function is written as 

where the argument of J in the integral arises from the invariant mass of the internal quark, (P2 — k + I) 2 w 
— 2(x — / + /P 1 + )Pi • Pi. Performing the integration over and It, we derive the evolution equation, 

dx 2n J x C~ x y ' V ' 

with the integration variable £ corresponding to the plus component Z + . 
Employing the variable change £ = xj z, Eq. ( |2^ ) becomes 

x dJ^)_ = J(0) + a. r 1 dz 

ax Zir J x z 

where the nonperturbative inhomogenious term j(°\x) absorbs the soft pole from £ — > x (z — > 1). After removing 
this pole, the integral in the second term is infrared finite. It is interesting to observe that Eq. ( p4| ) is similar to 
the evolution equation for unintegrated parton distribution functions involved in inclusive QCD processes ]20[| , which 
resums the same double logarithm a s In 2 x. 

To solve Eq. (p3|), we perform the Mellin transformation from the momentum fraction (x) space to the moment 
(JV) space: 

J(N)= [ dx{l- x^^Jix) . (25) 
Jo 

The small x region then corresponds to the large N region. The left-hand side of Eq. ( p3|) leads to 

1 dx(l - xf^x^-^- = -NJ(N) + (N- 1)J(N - 1) « --^7 \NJ(N)] , (26) 
dx dN L 

which is reasonable under the approximation dN sa AN = 1. The right-hand side is given by 

r l r l . /1 t \ N-l 



27T c F y Q ded-o-v(o y o — - ?^^) - (*) 

where we have exchanged the sequence of the integrations over x and over £, made the variable change x — £w, 
and kept only the jV-dependent term. The trick for extracting the N dependence from the integral over w in 
Eq. (pTl) is to consider its derivative with respect to N, and to insert the large N approximation for the <5-function, 
limjv^oo JV(l - y)"- 1 = S(y) with y = £(1 - tu)/(l - fte). 

Equating Eqs. ( p6|) and (|2^), we have the evolution equation in the moment space, 

N ^!m L = ~ J{N) ~ ^ c ^ Nj ( N ) • ( 28 ) 



whose solution is the Sudakov factor, 



1 / 1 , 2 



J(iV) = -expl -- 7jr ]n'J\M , (29) 

with the anomalous dimension 7^ = a s Cp/-K. The factor l/N can be regarded as the Mellin transformation of the 
initial condition j(°> = 1. Comparing Eq. ( p9[ ) with the resummation for an inclusive jet PJlCj], the resummation for 
an exclusive jet is suppressed by a factor l/N, and the corresponding anomalous dimension is down by a factor 2. 
By means of the inverse Mellin transformation, the jet function in the momentum fraction space is written as 

rc+ioa ijy 

J(x)= (l-a)-^), (30) 

Jc—ioo ^ 7r * 



G 



with c being an arbitrary real constant larger than all the real parts of thepoles of the integrand. Since J(N) contains 
a branch cut along the negative real axis on the complex TV plane, Eq. (|30| ) reduces to 

J(x) = -exp(^a s C F ) J°° ^(l-xrPWsinQa s C F ^exp(-g^ 2 ) , (31) 

where the variable change N = exp(i + in) (N = exp(i — in)) has been adopted for the piece of contour above 
(below) the branch cut as shown in Fig. 3. Note that the above expression holds only for a s > 0. For a s — 0, the 
residue associated with the N = pole should be included. It is trivial to check that J(x) is normalized to unity, 
fj(x)dx = J(l) = l & 

Obviously, Eq. (O) vanishes at x — > 0, because the integrand is an odd function in t, and at x — > 1 due to the 
factor (1 — x) cxp ^'. The latter property is the consequence of the extrapolation of the Sudakov factor to the low N 
region. Moreover, Eq. ( |3l| ) provides suppression near the end point x — > 0, which is stronger than any power of x. 
This is understood from vanishing of all the derivatives of Eq. (J3l]) with respect to x at x — ^ 0. For example, the first 
derivative gives 

^^)U- = -exp(|a s C F+ ^-) £^sin(la s C.t)exp(-^C^) , (32) 

where the variable change t — > t + 2n/ (a s C F ) has been made. The integrand in Eq. (|32] ) is also an odd function in t, 
and the integral diminishes. To the accuracy of the next-to-leading logarithms, the running of the coupling constant 
a s should be taken into account, and Eq. (|3l]) will be modified. However, the above features remain. 



V. SEMILEPTONIC DECAY 



In this section we extend the above formalism to the semileptonic decay B — > n I v in the fast recoil region of the 
pion. The B meson momentum P\ is the same as in the decay B — > "/W, and the pion momentum Pi is the same 
as the photon momentum. Leading- twist factorization theorem for the B — > n form factor f + (q 2 ) in the standard 
definition has been proved in Q] , 

/+(<? 2 ) = ^2 <f>m(xi) ® H m (x 1 ,x 2 ,r]) <g> (f> n (x 2 ) , (33) 

171= + 

which holds in the region with x\ ~ 0(K/Mb) and with X2 ~ O(l). 

Since Fig. 4(a), proportional to l/(a;ix|), is more singular at small x 2 , we consider the end-point region with 
X2 ~ 0(A/Mb), where the internal b quark propagator scales like l/(AMs). Part of 0(a s ) corrections to Fig. 4(a) 
are shown in Figs. 4(c)-(e), among which Fig. 4(c) generates the double logarithm a s In 2 X2 from the collinear region 
with the loop momentum parallel to P2. Figures 4(d) and 4(e) involve only the single logarithm. We emphasize that 
the double logarithm discussed here for Fig. 4(a) differs from that in jj], which concerns a s In 2 x%. As stated above, 
Fig. 4(a) is less singular at x\ — > 0, and a s In 2 X\ is not relevant. The double logarithm in Fig. 4(c) is simply extracted 
by eikonalizing the light quark line, which flows into the pion, and the b quark line according to Eq. (^|) : 



II J (2tt) 4 n+ • (/ - x 2 P 2 )n- ■ ll 2 ' 

= -^C F \n 2 x 2 + -- - . (34) 

Similarly, the collinear pole has been regularized by including a small plus component for n_ . It is observed that the 
double logarithm in Eq. (|34|) is the same as in Eq. (Q). 

The all-order factorization of the jet function from the decay B — > jlu can be proved following the procedures 
in Sec. III. It is also easy to show that this jet function obeys the evolution equation in Eq. ( f23[ ) by identifying the 
correspondence of £2^2 and n + to k and P2 in Sec. IV, respectively. Hence, the threshold resummation leads to a 
result the same as Eq. (|3l|). That is, the Sudakov factor is universal. 

The analysis for Fig. 4(b) is similar to that for the decay B — > jlu. In the end-point region with x\ ~ 0(A 2 /M|), 
additional collinear divergences associated with the internal light quark are produced in Figs. 4(f)-4(h). Figure 
4(f) gives the double logarithm a s ln 2 a;i, whose factorization is the same as of Fig. 1(c). Note that Fig. 5 has 
been considered as the dominant source of the double logarithm a s In 2 x\ associated with Fig. 4(b) Q). It will be 
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demonstrated explicitly in the Appendix A that Fig. 5 in fact contains only the single logarithm a s \nxi. The result 
in PJ might be attributed to an inappropriate collinear approximation. Therefore, the universality of the double- 
logarithm resummation indeed holds. 

At last, we investigate the resummation effect on the B — > ir form factor f + (q 2 ), whose factorization formula is 
given by 



/ + (? 2 ) 



dxidx2<t>B{xi) 



X x^ 



V 2 M 2 B N C 

with N c = 3 being the number of colors. We employ the models 

<t> n (x) = 6x(l - x){\ + 0.66[5(1 - 2a;) 2 - 1]} , 

(f>B (x) 



J(X2) n . s. J(xi) , ' 

(1 + x 2 r\) (1 - rj) 



X\X 2 



<t>A x ?) 



(35) 



1, 



b+(x) - 4>_{x)] = N B y/x{l-x) 



exp 



xM E 



(36) 



with the shape parameter cub and the normalization constant Nb determined by j 4>B{x)dx = 1. Another B meson 
distribution amplitude 4>b{x) — \4>+(x) — </>_(x)]/2 is power-suppressed, since its first nonvanishing moment starts 
with 0(X/Mb) p6[ . It has been confirmed, contrary to the conclusion in pi} |, that the contribution from (j>B is 
less important than that from (fis psfl . For the purpose of this work, it is sufficient to work on a single B meson 
distribution amplitude. If J(x) is excluded, the first term in Eq. ( |35"| ) is logarithmically divergent. With the threshold 
resummation, the form factor is calculable without introducing any infrared cutoffs 

Choosing a s = 0.4, M B = 5.28 GeV, and the decay constants f B = 190 MeV and / w = 130 MeV, we obtain f + (q 2 ) 
for uj b = 0.3(0.4) GeV, 



/+(0) = 0.15(0.12) , / + (2 GeV 2 ) = 0.18(0.13) , /+(4 GeV 2 ) = 0.21(0.15) 



(37) 



The variation of / + for a s — 0.3-0.5 is less than 10%, because the change of J is compensated by that of a s in 
the overall coefficient in Eq. (|3^). The difference between the above values and the expected one /+(0) ~ 0.3 can 
be resolved by taking into account higher-twsit contributions Jig] , Including two-parton twist-3 pion distribution 
amplitudes, which are finite at the end point |22]| , the singularities in the B — > ir form factor become linear. In this 
case threshold resummation is even more crucial. 



VI. CONCLUSION 



In this paper we have shown that the double logarithmic corrections a s In 2 x appear in exclusive B meson decays. 
When the end-point region with a momentum fraction x — > is important, these double logarithms need to be 
organized to all orders in order to justify perturbative expansion. The double logarithms, associated with internal 
particles which almost go on mass shell, are of the collinear nature. The PQCD factorization theorem for exclusive 
B meson decays at the end point then demands the introduction of a quark jet function into decay amplitudes. The 
factorization of the jet function from the B — > jW mode has been proved rigorously. The proof for the factorization 
from other modes, such as the semileptonic decay B — > nW, is similar. It has been shown that the jet function is 
defined as a matrix element of a quark field attached by a Wilson line, based on which threshold resummaiton can 
be performed. It is interesting to achieve the above proof in the framework of the soft-collinear effective theory |2(J . 

We have derived the evolution equation for the jet function, which is solvable in the large N limit in the Mellin space. 
The solution is a universal, i.e., process- independent, Sudakov factor. The qualitative behavior of this Sudakov factor 
has analyzed and found to decrease quickly at x — > 0. It has been demonstrated that Sudakov suppression is strong 
enough to smear the end-point singularity in the B — > ir form factor and leads to reasonable predictions. We conclude 
that in a self-consistent PQCD analysis of the heavy-to-light transition form factors, the end-point singularities do 
not exist. In a future work we shall extend the formalism developed here to more complicated nonleptonic B meson 
decays. Threshold resummation for various topologies of decay amplitudes, such as annihilation and nonfactorizable 
ones, will be studied. 

I thank M. Beneke, Y.Y. Keum, T. Kurimoto, D. Pirjol, A.I. Sanda, and G. Sterman for useful discussions. The 
work was supported in part by the National Science Council of R.O.C. under Grant No. NSC-90-2112-M-001-077, 
by the National Center for Theoretical Sciences of R.O.C, by Grant-in Aid for Special Project Research (Physics of 
CP violation), and by Grant-in Aid for Scientific Exchange from the Ministry of Education, Science and Culture of 
Japan. 



8 



APPENDIX A: TRIPLE-GLUON CORRECTION 



In this Appendix we show that Fig. 5 produces only the single logarithm a s lnxi. The loop integral is written as 



dH 



rtr 



75" 



75 P2 /P 2 - fa+ / 



(2tt)< 
Pi- M+ I 



4A C 
M B 



-IP 



(p 2 -fc 1 + 27 



7a- 



pa/36 jabc^^rpa'j^brj^c'^ 

(h - x 2 P 2 - l) 2 l 2 (x 2 P 2 - fcl) 2 



with the triple-gluon vertex, 



(2x 2 P 2 - 2Ai + i)V 4 + (^i - z 2 P 2 - 2l)°g ap + x 2 P 2 + ki) p g- 



8a@ 



\Pa8 



(Al) 



(A2) 



the anti-symmetric tensor 



fabc 



and the color matrices T a ' b ' c . 



We concentrate on the frist term in Eq. (A2), since the calculation of the seco nd and third terms is similar. If the 
loop momentum I in the quark propogator 1/(P 2 — k\ +1) is neg lected, Eq. @ gives the double logarithm a s In X\ 
as obtained in [^j. However, this approximation is not appropriate for the collection of the collinear divergence, for 
which I is of the same order as P 2 . Keeping the I dependence, and using the identity f abc t r (T a T h T c ) — 6i, we have 



h 



ig 4 tr{/P 2 frg^ Pi) 
4A C (x 2 P 2 -fc 1 ) 2 

dH v ■ {2x 2 P 2 - 2fci + I) 

(2tt) 4 v ■ l(P 2 -k!+ l) 2 (h - x 2 P 2 - l) 2 l 2 



(A3) 



The contribution proportional to M B in Eq. ( Al) is suppressed by a power of x\ compared to Eq. (|A3|). For a similar 
reason, the term v ■ k% in Eq. (A3) is also negligible. The term v ■ I in the numerator does not lead to a collinear 
divergence, because the poles of the corresponding denominator (P 2 — k\ + l) 2 (ki — x 2 P 2 — l) 2 l 2 in the l~ plane are 
not of the pinched type. 

Doing the contour integral with the residue l~ = —l + — ie (for l + < 0), Eq. (|A3|) becomes 



h - H 4b 



(2tt) 3 C f 



dP 



d 2 h 



x 2 p 2 {p 2 -K 



{l 2 T + 2l+ 2 



2P 2 ~l+ 



where 



H 



(0) M 
4b 



2P 2 -fc+)(Z 2 + 21+ 2 - 2x 2 P 2 l+ + 2x 2 P 2 k+)(l 2 r + 21+ 2 ) 
g 2 C F tr(JP 2 Ml" A) 



SN c (x 2 P 2 -k 1 ) 2 (P 2 -k 1 ) 2 ' 
is the lowest-order amplitude from Fig. 4(b). The integration over It leads to 

f]l + /M E 



(A4) 



(A5) 



/> OO 

HiT^c A x 2 k+ y o di- 



ln(/+ 2 /M 2 



xirj) 



(I - x 2 )(l+ + k+) 2 



ln(l+ 2 /M 2 + x 2 r,l+/M B + Xl x 2V ) \n(l+ 2 /M 2 ) 



x 2 (l - x 2 )(l+ + k?) 



+ \2 



x 2 {l+ + k+f 



(A6) 



with the color factor Ca = 3. 

The terms l +2 /M B in the first two logarithms, givi ng contributions suppressed by a power of X\ compared to those 
from l + /Mb, can be dropped. Each term in Eq. (A6) then produces only the single logarithms a s hix\: 



J i - H 4b 4^ CaX2 



lnxi 



hixi 

l-x 2 x 2 (l-x 2 ) 



x 2 



H lb -^C A \riXi 



(A7) 







In the above expression the cons tant terms independent of hxxx are not displayed. It can be shown in a similar way 
that the other two terms in Eq. (A2) involve only the single logarithms. 
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Figure Captions 

Fig. 1. (a) and (b) Lowest-order diagrams for the decay B — > jlD. (c)-(h) 0(a s ) corrections to Fig. 1(a). 

Fig. 2. (a) The Ward identity, (b) Factorization of the 0{a s ) diagrams as a consequence of (a), (c) Factorization of 
the B — ► decay amplitude in the end-point region. 

Fig. 3. Contour for the inverse Mellin transformation of the jet function. 

Fig. 4. (a) and (b) Lowest-order diagrams for the decay B — > nlD. (c)-(e) [(f)-(h)] Part of 0(a s ) corrections to 
Fig. 4(a) [Fig. 4(b)]. 

Fig. 5. Triple-gluon correction to Fig. 4(b). 
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